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Abstract 

We calculate the masses of the lowest lying eigenstates of improved SU(2), SU(3), SU(4) and 
SU(5) Hamiltonian lattice gauge theory (LGT) in 2+1 dimensions using an analytic variational 
approach. The ground state is approximated by a one plaquette trial state and mass gaps are 
calculated in the symmetric and antisymmetric sectors by minimising over a suitable basis of 
rectangular states. Analytic techniques are developed to handle the group integrals arising in the 
calculation. 
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I. INTRODUCTION 



In this paper we calculate the lowest lying glueball masses for SU(iV) LGT in 2+1 di- 
mensions, with N = 2, 3, 4 and 5, extending an earlier paper which considered N = 2 
and 3 fl. We use Kogut-Susskind classically improved and tadpole improved Hamil- 
tonians |3j in their calculation and develop analytic techniques for the calculation of the 
required expectation values. 

The outline of this paper is as follows. In Section |H] we develop analytic techniques for 
use in the calculation of certain group integrals for general SU(iV). When used as generating 
functions these integrals allow an analytic treatment of the matrix elements appearing in 
later sections. After defining our notation in Section UTT1 we fix the variational vacuum wave 
function in Section HVl Following that we calculate lattice specific heats in Section IVl before 
studying SU(iV) glueball masses in Section IVTl Section IVIIII contains our conclusions and a 
discussion of further work. 

II. ANALYTIC TECHNIQUES FOR SU(iV) 
A. The special cases of SU(2) and SU(3) 

In this section we recall two results which are useful in the calculation of group integrals 
for the special cases of SU(2) and SU(3). For further details we refer the reader to Ref. Q. 
Before presenting these results the concept of group integration needs to be introduced. For 
this purpose we introduce the one plaquette trial state, which we will use to simulate the 
ground, or (perturbed) vacuum, state with energy E Q , 



Here, |0) is the strong coupling vacuum defined by £?(x)\0) = for all i, x and a = 
1,2,..., A^ 2 — 1. £f{x) is the lattice chromoelectric field on the directed link running from 
the lattice site labelled by x to the site labelled by x + ai. The directed square denotes the 
traced ordered product of link operators, Ui(x), around an elementary square, or plaquette, 
of the lattice, 




(1) 




(2) 
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where a is the lattice spacing. With this notation understood, we can write the expectation 
value of a plaquette as an SU(iV) group integral as follows 



Q) : = 



(3) 



Here the products over I and extend over all links on the lattice, while the products over p 
and p' extend over all plaquettes on the lattice. In Eq. (jSJ) we introduced the notation Z\(p) 
to denote the trace of plaquette p. For each integral in Eq. Q the integration measure is 



given by the 



measure) 



3aar measure (also called the invariant measure and less commonly the Hurwitz 



For any compact group G, the Haar measure is the unique measure dll on 



G which is left and right invariant, 

j dUf{U) = [ dUf(VU) = [ dUf{UV) MV E G (4) 
Jg Jg Jg 

and normalised, 

dU = 1. (5) 



G 

In Eq. (j3J) / is an arbitrary function over G. 

In 2+1 dimensions the variables in Eq. Q can be changed from links to plaquettes 
with unit Jacobian p. The plaquettes then become independent variables allowing the 
cancellation of all but one group integral in Eq. J3J). All that remains is 

Here, V is a plaquette variable with TrV = -Zi(O). For the case of SU(2), analytic expressions 
for the plaquette expectation value in terms of modified Bessel functions have been used in 
variational calculations for almost 20 years. A key result []| is 

dUe c1rU = l h (2 C ). (7) 

SU(2) c 

Here I n is the n-th order modified Bessel function of the first kind defined, for all integers 
7i, by 

°° x 2k+n 

k=0 v 1 



In an earlier paper, we showed that the corresponding SU(3) result follows simply from 
a paper of Eriksson, Svartholm and Skagerstam |8j, with the result: 

In the next section we derive the general SU(iV) result as well as a selection of other 
useful SU(iV) group integrals. 



B. The generalisation to SU(iV) 

1 . Introduction 

Much work has been carried out on the topic of integration over the classical compact 
groups. The subject has been studied in great depth in the context of random matrices and 
combinatorics. Many analytic results in terms of determinants are available for integrals of 
various functions over unitary, orthogonal and symplectic groups [9]. Unfortunately similar 
results for SU(iV) are not to our knowledge available. Their primary use has been in the 
study of Ulam's problem concerning the distribution of the length of the longest increasing 
subsequence in permutation groups Q, till- Connections between random permutations 
and Young tableaux allow an interesting approach to combinatorial problems involv- 
ing Young tableaux. A problem of particular interest is the counting of Young tableaux of 
bounded height which is closely related to the problem of counting singlets in prod- 
uct representations mentioned in Section III Al Group integrals similar to those needed in 
this paper have also appeared in studies of the distributions of the eigenvalues of random 
matrices jlj, Q • 

In the context of LGT not much has work been done in the last 20 years on the subject 
of group integration. The last significant development was due to Creutz who developed 



16] 



using link variables 



a diagrammatic technique for calculating specific SU(iV) integrals 
This technique allows strong coupling matrix elements to be calculated for SU(iV) \7] and 
has more recently been used in the loop formulation of quantum gravity where spin networks 
are of interest [13, U, boj ■ 

In Sections III B 21 and III B 31 we extend the results of Section III Al to calculate two im- 
portant SU(A^) integrals. As generating functions these integrals allow the evaluation of all 



expectation values appearing in variational calculations of SU(iV) glueball masses in 2+1 
dimensions. To calculate these generating functions we work with plaquette variables and 
make use of techniques which have become standard practice in the fields of random matri- 
ces and combinatorics. In Section 111 B 21 we derive a generating function which allows the 
calculation of integrals of the form 

/ dU(TiU) m JwTe c ^ u+TlUt \ (10) 

JsV(N) 

The work in Section 111 B 31 generalises the generating function of Section 111 B 21 allowing the 
calculation of more complicated integrals of the form 

dU [Tr(^)] m e c ( Trf/+TWt ). (11) 

SU(JV) 

For each integral considered the approach is the same and proceeds as follows. We start 
with a calculation of a U(iV) integral. For example in Section Til B 21 we calculate 

G v{N) (c,d) = [ dUe^ u+dT * u \ (12) 

JU(N) 

This is a generalisation of G\j(n)(c,c), an integral first calculated by Kogut, Snow and 



Stone |2l(. We then make use of a result of Brower, Rossi and Tan [22j to extend the 
U(iV) integral to SU(iV) by building the restriction, det U = 1 for all U G SU(iV), into the 
integration measure. In this way SU(iV) generating functions can be obtained as sums of 
determinants whose entries are modified Bessel functions of the first kind. 



2. A simple integral 



In this section we introduce a useful technique for performing SU(iV) integrals. We start 
with the U(iV) integral of Eq. (|12j) and calculate it using a technique which has become 
standard practice in the study of random matrices and combinatorics. 

Since the Haar measure is left and right invariant (see Eq. (0J)) we can diagonalise U 
inside the integral as 



/ 



,i<j>a 



u = v 











\ 



(13) 



In terms of the set of variables {<fik}k=i the Haar measure factors as dU = d/j,((j))dV 



Since the integrand is independent of V, the V integral can be carried out trivially using 
the normalisation of the Haar measure given by Eq. 
Making use of the Weyl parameterisation for U(iV) 

= n|W)i 2 , (14) 



00, 



where A(0) is the Vandermonde determinant, with implicit sums over repeated indices 
understood, 



A(0) 



1 



icf>l(N-h) i4>2(N-i 2 ) . . . i<f} N (N-i N ) 



(15) 



we can express the U(iV) generating function as follows 



G\j(n)(c, d) 



2~ 



2tt 



2tt 



27T 



exp 



N 



^(ce** + de 



i=l 



|A(0)p 



(16) 



In Eq. (fl5|) . e^...^ is the totally antisymmetric Levi-Civita tensor defined to be 1 if 
{ii,..., i n } is an even permutation of {1, 2, . . . , n}, —1 if it is an odd permutation and 
otherwise (i.e. if an index is repeated). Substituting Eq. (fTKJ) in Eq. (fTH|) gives, 

N ,2n 



G\J(N)( c ,d) — j^^i 1 i 2 ...i N £jij2—jN 



TT / * ^ exp [z(j fe - * fc )0* + ce^ + de~^] . (17) 



To simplify this further we need an expression for the integral, 

" 2n dx 



9n(c, d) 



2tt 



exp(mx + ce lx + de 



(18) 



which is easily handled by expanding the integrand in Taylor series in c and d, 

OO OO J. 



c d f dx „• 



k=0 1=0 



ix(k— l+n) 



k\(k + n)\ 

I n (2Vcd) . 



k=0 
C 



(19) 



Making use of Eq. (fTT3|) in Eq. (fTTj) gives an expression for Gxjrm(c, d) as a Toeplitz deter 
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minant (a determinant of a matrix whose (i, ?)-th entry depends only on j — i), 

N 



Gv(N)(c, d) 



«JV c 'JlJ2.--JJV 



1 

AH 

1 



k=l 

k=l 



2\^tj 



det 



(20) 

l<i,3<N 

Here the quantities inside the determinant are to be interpreted as the (i,j)-th entry of an 
N x N matrix. Now to calculate the corresponding SU(iV) result the restriction det {7 = 1, 
which is equivalent to J2 k=1 <\>k = mod 2ti in terms of the (fik variables, must be built into 
the integration measure. To do this we follow Brower, Rossi and Tan |22j and incorporate 
the following delta function in the integrand of Eq. (JTfij) : 

(N \ oo / N \ 

^0 fc -Omod27r] = ^ 2ix5 I ^ <p k - 2nm J . (21) 
fc=l / m=— oo \fc=l / 

This is most conveniently introduced into the integral via its Fourier transform, 

oo / TV 

exp I zm ^ <^fc 

m=— oo \ fc=l 



(22) 



To obtain the SU(iV) integral from the corresponding U(iV) result the modification is there- 
fore trivial. Including Eq. (j22j) in the integrand of Eq. (fT^j) leads to the general SU(iV) 
result, 



G 



SU(iV) 



'c, d) 



dUe 



fTrU+dTrW 



SU(iV) 



^ det [9m+j-i{c, d))^ 



<N ■ 



(23) 



This expression can be manipulated to factor the d/c dependence out of the determinant as 
follows, 



det 



(l+j-i)/2 



h+j-i (2Vcdj 



l<i,j<N 



1 



c) 



N 

X \\ h+Om-i 
m=l 

-j det 



2Vcd 



2Vcd 



l<i,j<N 



(24) 



Making use of this result in Eq. leads to the SU(iV) generating function 



°° fd\ lN/2 
G S v(N){c,d) = ( _ ) det 

l=-oo ^ ' 



(25) 



For the case of SU(2) we can show that this reduces to the standard result of Arisue given 
by Eq. (J7J). To do this we need the recurrence relation for modified Bessel functions of the 
first kind, 



In-l{x) - In+l{x) = —I n (x). 

X 



(26) 



Recall that for SU(2) the Mandelstam constraint is Tr£7 = TrC/t, so the case Gsu(2)(c, c) can 
be considered without loss of generality; 



G SU (2) (c,c) = J2 [W 2 - Ii-i(2c)I l+1 (2c 



l=—oo 



J (4c)-/ 2 (4c). 



(27) 



Here we have used the standard addition formula 
ploying the recurrence relation of Eq. (|26[) gives 

1 



24j for modified Bessel functions. Em- 



Csu(2)(c, c) 



2c 



h (4c), 



(28) 



which is the standard result of Arisue given by Eq. (|7|). 

With an analytic form for SU(iV) in hand we can attempt to find simpler expressions for 
Gsu(N) (c, d) analogous to Eq. (|2*8*j) . To our knowledge no general formulas are available for 
the simplification of the determinants appearing in Eq. (|25j) . Without such formulas we can 
resort to the crude method of analysing series expansions and comparing them with known 
expansions of closed form expressions. Since the determinants appearing in the generating 
function are nothing more than products of modified Bessel functions we expect that if a 
closed form expression for the general SU(iV) generating function exists, it will involve the 
generalised hypergeometric function. With this approach we have limited success. The SU(3) 
result of Eq. (jHJ) is recovered numerically but the SU(4) result does simplify analytically. 

When analysing the series expansion of Gsu(4) ( c ? c ) we notice that it takes the form of a 
generalised hypergeometric function. In particular we find the following result: 

Gsu(4) (c,c) =2^3 " ' :16<- 



3,4,5 



(29) 
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Here the generalised hypergeometric function is defined by 



p r q 



Oi, Q2, ■ ■ ■ , «% 



61,62 



X 



,b q 



k=0 



(ai)fc(Q2)fc • • • {a P )k X k 
{bi)k{b 2 )k ■ ■ ■ (b q )k k\ 



(30) 



where (x)k = x(x + 1) . . . (x + k — 1) is the rising factorial or Pochhammer symbol. In 
addition to Eq. (|29|) we find the following results for matrix elements derived from the 
SU(4) generating function: 



C2-P3 



(Zi)- 



5 1 

2 ' 2 ;16c 2 
4,5,6 



2-^3 



3 5 
2' 2 

3,4,5 



16c 2 



(31) 



3„2 vp 
9 C 2-^3 



5 7 




7 9 




9 11 


2 ' 2 ;16c 2 


+ |c 4 2 F 3 


2 ' 2 ;16c 2 


+ ^ 2 F 3 


2 ' 2 ; 16c 2 


5,6,7 




_6,7,8 


7,8,9 



2-^3 



3 5 
2' 2 

3,4,5 



16c 2 



and 



(ZiZi)-- 



2F3 


3 5 

2 ' 2 ;16c 2 
3,5,6 


+ |c 2 2 F 3 


5 7 

2 ' 2 ;16c 2 
_4,6,7 


+ |c 4 2F3 


7 9 

2 ' 2 ;16c 2 
_5,7,8 


2F3 


3 5 

2 ' 2 ;16c 2 
3,4,5 





(32) 



(33) 



We stress that these results are nothing more than observations based on series expansions. 
Despite some effort analogous expressions for N > 4 have not been found. 

The generating functions, Gs\j(n)(c, d) and Gu(iv)( c 5 d), are not only of interest in Hamil- 
tonian LGT. By differentiating Eq. (|23|) appropriately with respect to c and d and afterwards 
setting c and d to zero, we obtain the number of singlets in a given product representation 
of SU(iV). This was discussed for the special case of SU(3) in Section III Al We now consider 
the general case in the calculation of Tfc(n); the number of singlets in the SU(fc) product 
representation, 



(k <g> k) <g> • • • ® (k <g> k) . 



(34) 



n 



As a group integral Tk(n) is given by 

T k (n) = 



dU{\TrU\ 



2\n 



(35) 



SU(fc) 



Integrals of this kind are studied in combinatorics, in particular the study of increasing 
subsequences of permutations. An increasing subsequence is a sequence i 1 < i 2 < ■ ■ ■ < i m 
such that 7r(«i) < 7r(z 2 ) < ■ • • < 7t(i m ), where 7r is a permutation of {1, 2, . . . , k}. It has 
been shown that the number of permutations 7r of {1, 2, . . . , k} such that 7r has no increasing 
subsequence of length greater than n is T k {n) [10|. In addition it is possible to prove that 
Tfc(ra) is the number of pairs of Young tableaux of size k and maximum height n via the 
Schensted correspondence |ld 125]. 

Making use of Eq. (J25|) we see that only the I = term contributes to Tk(n). Letting 
x = cd we have: 



T k (n) = (d c d d ) n det Ij^Vcd) 
. d 



nl 



dx r 



J l<i,j<k 

det [1^2^)}^.. 



c=d=0 



j<k 



(36) 



x=0 



Hence the generating function for T k {n) is given by 

E^l^ = det[J^(2V5)] 



n=0 



nl 



l<i,j<k 



(37) 



a result first deduced by Gessel |l3(. The first few Tk{n) sequences are available as A072131, 
A072132, A072133 and A072167 in Sloane's on-line encyclopedia of integer sequences 



3. A more complicated integral 



We now move on to the more complicated integral 



H m (c,d) = [ dUexp [c(TrU + TiU j ) + dTr(U m )] Vm 6 

JSV(N) 



(38) 



This integral is of interest as a generating function for the calculation of integrals such as 



/ dUTr(U m )e 

JSV(N) 



m^cCIW+TK/t) 



(39) 



For the simple case of SU(3) we can use the Mandelstam constraint, Ti(U 2 ) = (Till) 2 — 
2TtW , to reduce such integrals to those obtainable from Gsu(at)(c, d) . However for higher 
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dimensional gauge groups not all trace variables can be written in terms of Till and TrW . 
For these gauge groups one must introduce the generating function, H m (c,d), to calculate 
integrals similar to Eq. (|39p. 

To calculate H m (c, d) we start with the corresponding U(iV) generating function h m (c, d) 
and follow the procedure of Section 111 B 21 to obtain 



h m {c,d) = [ dUexp [c(TiU + TiU r ) + dTi(U m )] 

JV(N) 

h ... jN J] / ^ ex P [*& - ^)<Pk + 2ccos0 fc + de mi ^\ . (40) 



1 

M 



To proceed we need the following integral, 

f 27T dx °° b k f 27T dx 

/ — exptinx + a cos x + 6e im:c ) = — / — e 
./o 2tt ^ fc! y 2tt 



i(n+mk)x+a cos x 



fc=0 

Making use of Eq. (j41j) in Eq. ()4Uj) leads to the following expression for the U(iV) generating 
function, 

h m (c,d) = det [X m;j _i(c, d))^^ , (42) 

with 

°° d k 

A m;n (c, d) = ^ — J n+TOfc (2c). (43) 

fc=0 

Extending to SU(iV) following the prescription of Section III B 21 we arrive at the correspond- 
ing SU(iV) generating function, 

oo 

H rn (c,d) = ^ det [X m .i +j _i(c, d))^ ^ N . (44) 

l=— oo 

An example of an SU(iV) integral derived from this generating function is the following: 
dUTr{U m )e c 



jm \ c(TW+TWt) _ 9H m (c, d) 



su(;v) 9d 

r-, OO 

a 



d=0 

• (45) 

d=Q 



— ^ det [ J '+i-i( 2c ) + dli+j-i+m{2c)] 



Only two terms need to be kept in the fc-sum of Eq. (J43|) here because higher order powers 
of d vanish when the derivative with respect to d is taken and d set to zero. 
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III. ANALYTIC SU(iV) CALCULATIONS 



Preliminaries 



In this section we make use of the analytic results of Section |n] in variational calculations 
of glueball masses. The steps we take are as follows. We first fix the variational ground 
state by minimising the vacuum energy density. Having fixed the ground state we are then 
free to investigate the excited states. 

Before fixing the variational ground state we introduce the following convenient notation. 
We define the general order a 2 improved lattice Hamiltonian for pure SU(iV) gauge theory 
with coupling g 2 on a lattice with spacing a by 



H(k,u ) 



a 



as,* 

2N 
ag 2 



k)Si(x) + —S i (x)U i (x)S l (x + ai)Uj(x) 



Hp 



{(l + 4K)P„(aO - \ \R l3 {x) + R ji {x)\) 



(46) 



x,i<j 



where the plaquette and rectangle operators are given by 



Rij{x) 



(47) 



1 - -^ReP 
ulN U , 

The simplest lattice Hamiltonians derived in Ref . can be expressed in terms of 7i as follows. 



The Kogut-Susskind 



and 0(a 2 ) classically improved Hamiltonians [3] are given by H(0, 1 



and 71,(1/6, 1) respectively. The tadpole improved Hamiltonian |3( is given by TC(l/6,u ), 
where the mean link uq, which can be expressed in terms of the mean plaquette, 



UQ 



(48) 



is defined self-consistently as a function of (3 = N/g 2 as described in Section HV Al 
With this notation the vacuum energy density is given by 



a 



€0 = N p W 

(1-AC) 

2k/3 

~m4 



N 2 -l 

2(3 



2/3(1 + 4k) 



D 



(49) 



+ 2(1 + 3«)£, 
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where N p is the number of plaquettes on the lattice and the expectation values, as usual, 
are taken with respect to the one plaquette trial state defined in Eq. (0). The variational 
parameter, c, is fixed as a function of (3 by minimising the vacuum energy density. For the 
calculation of the expectation values we use the generating functions of Section |H] with all 
infinite sums truncated. Once the variational parameter is fixed the trial state is completely 
defined as a function of [3. 

IV. FIXING THE VARIATIONAL TRIAL STATE 

A. Introduction 

In this section we fix the SU(iV) trial state for 2 < iV < 5, making use of the generating 
functions derived in Section |HJ These generating functions allow the analytic calculation of 
the plaquette and rectangle expectation values appearing Eq. (J49|) . The approach we take 
is as follows. For the Kogut-Susskind and classically improved cases, we simply minimise eo 
for a given value of (3. The value of c at which eo takes its minimum defines c as a function 
of j3. The tadpole improved case is more complicated because the mean plaquette depends 
on the variational state, which is determined by minimising the energy density. The energy 
density however, depends on the mean plaquette. Such interdependence suggests the use 
of an iterative procedure for the calculation of the tadpole improved energy density. The 
approach we adopt is as follows. For a given j3 and starting value of Uq we minimise the 
energy density of Eq. (|49J) to fix the variational state |0o)- We then calculate a new mean 
plaquette value using this trial state and substitute in Eq. (l4*9"j) to obtain a new expression for 
the energy density. This process is iterated until convergence is achieved, typically between 
five and ten iterations. 

B. Results 

The results of the Kogut-Susskind, classically improved and tadpole improved SU(2), 
SU(3), SU(4) and SU(5) vacuum energy density calculations are shown in Fig. [TJ The 
corresponding variational parameters c(/3) are shown in Fig. 121 For SU(3) the generating 
function of Eq. (jHJ) is used to calculate the required plaquette and rectangle expectation 
values. The generating function of Eq. (J23j) is used for SU(4) and SU(5). 
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FIG. 1: Analytic calculation of the 2+1 dimensional unimproved, improved and tadpole improved 
vacuum energy density in units of l/(aN p ) for SU(2), SU(3), SU(4) and SU(5). 

The familiar strong and weak coupling behavior from variational calculations is observed 
in each case. The differing gradients for the improved and Kogut-Susskind c((3) for a given 
TV in the weak coupling limit highlight the fact that when using an improved Hamiltonian 
one is using a different renormalisation scheme to the unimproved case. 
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(c) SU(4) (d) SU(5) 

FIG. 2: Analytic calculation of the unimproved, improved and tadpole improved variational pa- 
rameter in 2+1 dimensions for SU(2), SU(3), SU(4) and SU(5). 

C. Dependence on truncation 

In practice the /c-sum appearing in the SU(3) generating function is truncated at a maxi- 
mum value k max . The dependence of the variational parameter on various truncations of the 
/c-sum is shown in Fig. |3J We see that convergence is achieved up to (3 ~ 13 when keeping 
20 terms. Further calculations show that when keeping 50 terms convergence up to f3 ~ 30 
is achieved. 

The Z-sum appearing in the general SU(iV) generating function of Eq. (J23j) is also trun- 
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FIG. 3: Analytic calculation of the unimproved SU(3) variational parameter in 2+1 dimensions, 
truncating the /c-sum of y(c, d) at k = k max . 
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FIG. 4: Analytic calculation of the unimproved SU(3) and SU(4) variational parameters in 2+1 
dimensions, truncating the /-sum of Eq. (|25|) at I = ±l max . 



cated in practice. We replace the infinite sum over I by a sum from — Z max to Z max . The 
dependence of the SU(3) and SU(4) variational parameters on Z max is shown in Fig. HJ From 
the graphs we see that convergence is achieved quickly as Z max increases for both SU(3) and 
SU(4). The results for Z max > 8 are barely distinguishable up to (3 = 80 with the scale used 
in the plots. 
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V. LATTICE SPECIFIC HEAT 



In addition to the vacuum energy density we can also calculate the lattice specific heat 



C V = -f*. (50, 



The results for SU(2), SU(3), SU(4) and SU(5) are shown in Fig. The SU(2) and SU(3) 
results are calculated with the aid of Eqs. (jZJ) and (|SJ) with the fc-sum of Eq. © truncated at 
fc mm = 50. The SU(4) and SU(5) results are obtained using Eq. (|25[) to calculate the required 
matrix elements. For these cases the infinite /-sum is truncated at Z max = 4, for which the 
generating function has converged on the range of couplings used. We recall that the location 
of the peak indicates the region of transition from strong to weak coupling j^j]. For an 
improved calculation one would expect the peak to be located at a smaller (3 (corresponding 
to a larger coupling) than for the equivalent unimproved calculation. We see that this is 
indeed the case for each example, with the tadpole improved Hamiltonian demonstrating 
the largest degree of improvement. 

VI. MASS GAPS 
A. Introduction 

Having fixed the one-plaquette vacuum wave function, in this section we turn to inves- 
tigating excited states. Our aim is to calculate the lowest lying energy eigenstates of the 
Hamiltonians described by Eq. (J46|) for SU(iV) with 2 < iV < 5. 

We follow Arisue (28 1 and expand the excited state |0i) in the basis consisting of all 
rectangular Wilson loops {\n, m)} n ^ 1 = {\1)}[™{- X that fit in a given square whose side 
length L max defines the order of the calculation. Enumerating the possible overlaps between 
rectangular loops is relatively simple and so a basis consisting of rectangular loops is an 
ideal starting point. However, for an accurate picture of the glueball spectrum we will need 
to extend the rectangular basis to include additional smaller area loops. Without such small 
area nonrect angular loops, it is possible that some of the lowest mass states will not appear 
in the variational calculation described here. In order to ensure the orthogonality of |</>o) 
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FIG. 5: The unimproved, improved and tadpole improved specific heat in 2+1 dimensions for 
SU(2), SU(3), SU(4) and SU(5). 



and we parameterise the excited state as follows 



i0i) = s 'i n ' m ) = Si ^> 



(51) 



n,m=l 



1=1 



with, 



|0 = ^[W^)-(^)>]l0o>. 



(52) 
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Here (Wi(x)) is the expectation value of Wi(x) with respect to the ground state \4>o) and 
the convenient label I — (n — l)L max + m has been defined to label the n x m rectangular 
state, \n, m). We define the particular form of W\(x) to reflect the symmetry of the sector 
we wish to consider. For SU(iV) we take Wi(x) = Tr[wi(x)±w](x)] for the symmetric (0 ++ ) 
and antisymmetric (0 ) sectors. To avoid over-decorated equations, the particular Wi(x) 
in use is to be deduced from the context. Here wi(x) is the rectangular Wilson loop joining 
the lattice sites x, x + nai, x + nai + maj and x + naj, with 



n 



I - 1 

L rn.n.'r. 



1 and m = I — L r 



I - 1 



(53) 



Here [k] denotes the integer part of k. In order to calculate excited state energies we minimise 
the mass gap (the difference between the excited state and ground state energies') over the 



basis defined by a particular order L r 
the matrices 

Niv 

where E is the ground state energy, and 



To do this we again follow Arisue 28 



{l\H-E \l'), 



and define 



(54) 



Div = -(iio = £[(^w(o))-wwr(^(o)) 



(55) 



Extending the calculation to the general improved Hamiltonian H(k, uq) following Ref. 



i 



gives 



Nw 



i.x x' 



+ " 



£?(aO,WiV) 



E?{x),w!(x') E?(x + ai),W v (0) 



[E?{x),W v {0)] 



(56) 



To minimise the mass gap over a basis of a given order we make use of following diagonali- 
sation technique (Q. We first diagonalise D, with 

S^DS = V 2 , (57) 

where V is diagonal. The n-th lowest eigenvalue of the modified Hamiltonian 



H' = VS^NSV, 



(58) 



then gives an estimate for the mass gap corresponding to the n-th lowest eigenvalue of the 
Hamiltonian, Am n . 
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B. Classification of states 



States constructed from only gluon degrees of freedom can be classified in the continuum 
by their J PC quantum numbers. The assignment of P and C quantum numbers is straight 
forward [30]. Care must be taken, however, when building states with particular continuum 
spins on the lattice. Difficulties arise when the continuous rotation group of the continuum 
is broken down to the group of lattice rotations. The most serious difficulty to arise is an 
ambiguity in the assignment of continuum spins to states built from lattice operators j^]. 
To give a specific example, suppose we construct a wave function, \(f>), on the lattice with 
lattice spin J = 0; a state built from Wilson loops which are unchanged by rotations of 
mi/ 2 for all integers n. As explained in Ref. I^, this state is not a pure J = state; it also 
contains J = 4, 8, . . . states. Using a variational approach we can obtain estimates of the 
lowest energy eigenvalues of states with lattice spin J = 0. When the continuum limit is 
taken we obtain estimates of the lowest continuum energy eigenvalues for the states with 
spin 0, 4, 8, . . .. 

In the Lagrangian approach it is possible to suppress the unwanted spin J ± 4, J ± 8, . . . 
states in a given spin J calculation. By "smearing" links, one can confidently construct 
lattice states which do not couple with the unwanted higher spin continuum states, at least 
for the lowest energy eigenvalues JscJ. The technique of "smearing" links has not to our 
knowledge been applied in the Hamiltonian approach. 

Another way to clarify ambiguities in spin assignment is to attempt to construct a J = 4 
state on the lattice, devoid of J = contributions, and similarly a J = state devoid of 
unwanted J = 4 contributions. The construction of an exact J = 4 state is impossible 
on a square lattice due to the unavailability of 7r/4 rotations. One can however attempt 
to construct states that are approximately symmetric under rotations by tt/4. In the La- 
grangian approach in 2+1 dimensions, it has been demonstrated, for the simple case of 
SU(2), that such states can be chosen on a square lattice and that the approximate it /4 
symmetry becomes exact in the continuum limit j^JJ. This technique is readily applicable 
in the Hamiltonian approach but has not yet been attempted. 

Thus in our Hamiltonian calculation, which incorporates only rectangular states, a lattice 
spin J state will correspond to a continuum state with spin J, J ± 4, J ± 8, . . .. Using a 
variational approach we can obtain estimates of the lowest mass states in the continuum 



20 



with these spin values. To improve the spin identification in the continuum beyond modulo 
4 requires more work. It will prove interesting to compare the masses calculated here to 
that of Teper 30j who has been careful to identify continuum spins correctly, at least for 
the lowest mass excitations. 



C. Calculating matrix elements 

Having described the minimisation process we now detail the calculation of the matrix 
elements and Du>. Our aim is to reduce Nu> and Du> to polynomials of one plaquette 
matrix elements. This, again has been done for the case of SU(2) by Arisue 28|. Here 
we retrace his calculations for the general case of SU(iV) and extend them to incorporate 
improved Hamiltonians. We start with D\\i. 

Taking elementary plaquettes as our independent variables, it is easy to show that the 
only non-zero contributions to Diy occur when the rectangles I and I 1 overlap. As an example 
of a contribution to Du>, consider ADu>] the case where Nmv plaquettes are contained by 
both rectangles (these are the overlap plaquettes) and Ni plaquettes are contained by the 
rectangle I. In order to calculate such matrix elements we rely heavily on the orthogonality 
properties of SU(iV) characters. We are interested in calculating SU(iV) integrals of the 
form 

J dU p e s ^ X r(U p V), (59) 

where U p is a SU(iV) plaquette variable and V is a product of any number of plaquettes not 
including U p . Here Xr{U) denotes the character corresponding to the representation r. For 
SU(2), r = 0, 1/2,3/2, . . . and for SU(3), r = (A,/i) where A denotes the number of boxes 
in the first row of the Young tableau describing the representation and fi is the number of 
boxes in the second row. Similarly, for SU(iV), r = (ri, . . . , rjv-i)- 
Performing a character expansion of the exponential in Eq. (joTJj) gives: 

J dU p e S ^ Xr {U p V) = Y.j d UpCr'Xr'(U p ) X r(U p V). (60) 

This is simply a generalisation of a Fourier expansion. Here, the coefficient cy is given by: 

= J dU pX AU p )e s ^\ (61) 
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Now, using the orthogonality relation, 

J dU P xAU p V) X r(U p ) = j<WXr(V0, (62) 

where d r is the dimension of the representation r, we obtain: 

J dU p e s ^ Xr (U p V) = jXriV) J dU pXr {U p )e s ^\ (63) 

This result allows us to integrate out a single plaquette from an extended Wilson loop in 
a given representation r. To complete the calculation we need to relate SU(iV) characters 
to the traces of group elements. This can be done using Weyl's character formula [s^. For 
SU(iV), according to Bars [3], the dimensions and characters corresponding to the first few 
representations are given by: 



Xi(U) 


= TrU 




di(U) 


= N 


X2(U) 


= \ [(Trt/) 2 


+ TrU 2 \ 


d 2 {U) 


= \n(n + i) 


Xn{U) 


= \ [(™) 2 


- Trf/ 2 ] 


dn(U) 


= In(n-i) 


X2l{U) 


= \ [(Trt/) 3 


- Trf/ 3 ] 


d 21 (U) 


= -(N-1)N(N+1 


Xi»-i(U) 






d lN -x(U) 


= N 



(64) 



Here we have adopted the convention of dropping all zeros in the character labels. The 
Mandelstam constraints for the gauge group in question allows all characters to be expressed 
in terms of a minimal set of trace variables. For example, for SU(3) we make use of the 
Mandelstam constraint, 

Tt(A 2 B) = TrATr(AB) - Tr^Trfi + Ti(A^B), (65) 

where A G SU(3) and B is any 3x3 matrix, to express all characters in terms of Tr£7 and 
Trf/t. For example, for the case of SU(3), Eq. (JMj) simplifies to 

Xi(U) = TrU di = 3 

Xn(U) = \ [{JiUf - Tr(f/ 2 )] = TrC/t d u = 3 

I (66) 
X2 (U) = - [(TrU) 2 + Tr(U 2 )} = (TrU) 2 - Trtf d 2 = 6 

Xai{U) = \ [(TiU) 3 - Ti(U 3 )] = TrC/TrC/t - 1 d 21 = 8 

However, for general SU(iV) the Mandelstam constraints are difficult to calculate. In what 
follows we will need expressions for TtU, TrW, Tr(U 2 ), (TiU) 2 , and Trf/Trf/t as linear 
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combinations of characters for SU(iV). It is possible to do this without the use of the 
Mandelstam constraint. Such expressions are necessary in order to make use of Eq. (JUHj) in 
the calculation of expectation values of trace variables. For TrU, TrU', Tr(U 2 ) and (TW) 2 
the necessary expressions are easily obtained by rearranging Eq. (J64)) . 



TrU 
(TrU) 2 
TrU 2 
Trtf 



Xi(U) 

X2(U)+xn(U) 
X2(U)- X ii(U) 
Xi»-i(U). 



(67) 



To express the remaining expression, Trf/TrfA, as a linear combination of characters is not 
as easily done. For SU(3) one can simply rearrange the expression for X2i{U) in Eq. (J66j) . 
For the general N case it is simplest to consider Young tableaux. In terms of Young tableaux 
we have 



TrUTrU ] = □ 



□ 
□ 

□ 



y N-i. 



(68) 



Performing the product representation decomposition gives 



□ 



□ 
□ 

□ 



\ N-l 



□ 
□ 

□ 



> N 



□ 
□ 



ViV-l 



(69) 



Converting back into the notation of characters and traces gives 



TrUTrU^ = 1+ X 2i»-*{ U )- 



(70) 



Eqn ffTOj) together with Eqs. (jBTj) and (jB^j) allow the analytic calculation of each contribu- 
tion to Dui f° r & h N. For the case of ADui described earlier, we have 

AD U , = ^F Zl (Nt + Ni> - 2N lnV ) [f^Nm) + F Zl2l (N ln i,)] - 4F Zl (Ni)F Zl (N v ), (71) 



23 



where the character integrals are defined by: 
FzAn) 



2 [N(N + 1] 



n-l 



(zl + z 2 >» + - 



Fz 2 (n) 



N(N + 1) 



n-l 



{zl + Z 2 ) n - 



N(N — 1) 
1 

N(N - 1) 



n-l 



n-l 



<^ 2 - z 3 >» 



[n)=l + 



~1 n-l 



(72) 
(73) 
(74) 
(75) 



_(N-1)(N + 1)_ 

Here we have made use of the notation, Z n := Tr {U n ), to denote the trace variables occupy- 
ing a single plaquette, U. The expectation values appearing in Eq. (fTHj) are easily expressed 
in terms of the generating functions of Eqs. (f23j) and (j44jl . Differentiating the generating 
functions appropriately gives 

1 dG S \j(N) 



(Zi) 
(Z 1 Z l ) 

(zl ± z 2 ) 



GsV(N) 


dc 


1 


d 2 G S \j(N) 


Csucat) 


dcdd 


1 


d 2 Gsv(N) 


Gsv(N) 


dc 2 



d=c 



d=c 



± 



1 dH 2 



d=c 



H 2 3d 



(76) 



d=0 



In practice, we do not need to calculate all of these matrix elements. We see from Eq. f!25j) 
that (Zf) and {ZiZ\) are related by 



(Zl) 



2 \ d\Z x ) 
2 dc 2 



(Z l Z 1 ). 



(77) 



This follows from the fact that a group integral does not depend on the choice of direction 
for the links. To be more precise, the result, 



dUf{U) = / dWf(W) 

SV(N) JSV(N) 



(78) 



SU(N) 



follows from the fact that dW and dU each define invariant Haar measures on SU(iV) which, 
by uniqueness, must be equal. 

We now move on to the calculation of Nu>. It is easy to show that the only non-zero 
contributions occur when there is at least one common link and an overlap between the 
rectangles. The improvement term (the second term in Eq. (jSHJ)) only contributes when 
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the two rectangles share at least two neighbouring links in a given direction. Consider the 
contribution ANu> to Nu> in which there are L\ common links and L2 common strings of 
two links in a given direction. Again we suppose N\ plaquettes are enclosed by rectangle I 
and that there are iVj n j/ common plaquettes. Making use of Eq. (J63j) and following Ref. Q 
we obtain 



AJV„/ = -F Zl (Ni+N v -2Nw) 



Fz 2 (AW) - \f z1 (N lnV )-N+±-F ZlZl (AW 



with 



(79) 



'1 - k)Li + —~L 



2- 



For the case of SU(3) we can simplify this using Eq. (jHHjl to 



AN IV = -FzM+Ny^Nw 



1 



-F Z 2 (N lnv ) - 2F Zl (N lnv ) - 3 + -F ZlZl (iV, 



ini' 



(80) 



(81) 



2 \ tl It / £J\ \ n It / 1 ^ 

Having determined individual contributions to Du> and Nu> , to complete their calculation 
the possible overlaps between states I and I' of a given type must be counted. 



D. Choosing an appropriate vacuum state 

In Section IIVI we calculated variational vacuum wave functions for pure SU(iV) gauge 
theory for iV = 2, 3, 4 and 5. Our motivation was to use these wave functions as inputs to 
calculations of SU(iV) massgaps. We obtained wave functions with a variational parameter 
that was proportional to (3 in the large j3 limit and 1 in the small (3 limit. However, this 
behaviour is incompatible with the exact continuum vacuum wave function. For a one- 
plaquette trial state, given by Eq. ((H), to be compatible with the exact continuum vacuum 
wave function in 2+1 dimensions, one must have c oc [3 2 in the large j3 limit jQ]. This result 
is independent of the dimension of the gauge group in question. For the case of SU(2) in 
2+1 dimensions, in the scaling region it has been shown that for compatibility with 
the exact SU(2) continuum vacuum wave function, we must use the Greensite vacuum wave 
function 



M = e*p ( E ReQ| |0>. 

V ^ x,i<j J 



(82) 



It is this vacuum wave function that was used by Arisue [28( in the calculation that we 
generalise in this paper. It would thus seem that using a variational wave function is not 
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appropriate in the calculation of massgaps. However, the scaling argument presented in 
Ref. |l| demonstrates that in a glueball mass scaling region the precise form of c(/3) is not 
important. For this reason we expect that the results presented here will agree with a 
calculation using the Greensite wave function. 

As it turns out, to calculate variational wave functions for large dimension gauge groups 
is cumbersome. Numerical precision becomes a factor in the minimisation of the energy 
density. This problem is magnified in the calculation of tadpole improved results. We thus 
abandon the use of variational wave functions in the calculation of massgaps beyond SU(5). 
Instead we make use of the one plaquette wave function of Eq. (JTJ) and define a simple 
dependence, c(/3), which most often will be c{(3) = (3. Calculations for N > 5 will be 
presented in a later publication. 

VII. SU(2), SU(3), SU(4) AND SU(5) MASSGAP RESULTS 

In this section we present glueball mass results for SU(iV) pure gauge theory in 2+1 
dimensions with N = 2, 3, 4 and 5. For each SU(3) calculation we keep 80 terms in the 
/c-sum of Eq. Q giving convergence up to (5 = 50. For N > 3 the truncation Z max = 20 
is used. The generation of Nu> and Dui and implementation of the minimisation process is 
accomplished with a Mathematica code. 

For the case of 2+1 dimensions we expect aAm/g 2 to become constant in the scaling 
region. The convergence of the massgaps with L max is illustrated in Fig. |H1 We notice that 
for N > 2 only small improvements to the scaling behaviour are gained by extending the 
calculation beyond order 8 on the range of couplings shown. This suggests that a more 
complicated basis (including, for example, nonrect angular loops) is required to simulate 
SU(iV) excited states with N > 2 than for the case of SU(2). 

In Fig. [7| results for the lowest lying glueball mass, calculated with Kogut-Susskind, 
improved and tadpole Hamiltonians, are shown. We see that aAmf/g 2 is approximated 
well by a constant, in very large scaling regions, for the lowest lying eigenstates for all N 
considered. The scaling behavior becomes significantly worse for the antisymmetric sector 
which is shown in Fig.|B]and for higher energy eigenvalues. This is because the simplistic form 
of our excited state wave function is not sufficient to reproduce the plaquette correlations 
required to simulate these higher order states. One would expect the simulation of higher 
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FIG. 6: The unimproved 2+1 dimensional symmetric massgaps for SU(2), SU(3), SU(4) and SU(5). 



order eigenstates to improve by including more complicated loops in our expansion basis or 
by using a more complicated ground state. The continuum limit excited states results for 
SU(iV) are given in Tables ITI-IVTTI 

For the unimproved SU(2) case, the masses of the lowest two eigenstates agree closely 
with the calculations of Arisue [2^ (respectively 2.056 ± 0.001 and 3.64 ± 0.03 in units of 
e 2 = g 2 / a) in which the Greensite vacuum wave function of Eq. ()82|) was used. This serves as 
a check on our counting in calculating the possible overlaps of excited states. Our calculation 
is in disagreement with that of Arisue at the third eigenstate, for which Arisue calculates a 
mass (5.15 ± 0.1)e 2 . Our fourth eigenstate is close in mass to Arisue's third and our third 
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eigenstate does not appear in his results. The reasons for this are not clear. 

The results for the SU(3) symmetric massgap (in units of e 2 ) are to be compared to calcu- 
lations by Luo and Chen 2. 15 ±0.06 34 1, Samuel 1.84 ±0.46 }a5 1 and Teper 2.40 ±0.02 Q- 
Our result of 3.26520 ± 0.00009 is considerably higher than all existing comparable results. 
By including more complicated loops in the expansion basis one would expect to reduce this 
estimate. This is emphasised by the fact that when using only square basis states our result 
is considerably higher. To explain the discrepancy between the our results and other's it 
is important to note that since we use a basis of rectangles we exclude the contribution of 
many non-rectangular small area diagrams that are included in the calculations of Teper and 
that of Luo and Chen. For this reason it may be the case that what we have interpretted as 
the lowest glueball mass in this paper may, in fact, be a higher order excited state. Teper 
has calculated the masses of the three lowest mass glueballs for SU(3) in the ++ sector [30]: 
0++, 0++* and 0++**, with the respective results, in units of e 2 : 2.40 ± 0.02, 3.606 ± 0.063 
and 4.55 ±0.11. It is interesting to note that our result is closer to Teper's first excited state. 
In the same study Teper also calculated glueball masses in the ++ sector for N = 4, 5 and 
6. The mass, in units of e 2 , of his ++ * state for SU(4) is 4.84 ± 0.12 and 5.99 ± 0.16 for 
SU(5). We notice that as iV is increased the results presented here move closer to the mass 
of Teper's ++ * state, with the improved results being closer than the unimproved. In fact 
for SU(5), the results presented here, improved and unimproved, are consistent with Teper's 
q++* mass _ This forces us to question the interpretation of the large (3 plateaux in Fig. 
as scaling regions for the lowest mass glueballs. It is possible that the minima present in 
Fig. in the small j3 region are possible scaling regions. It is possible that our vacuum wave 
function and minimisation basis are insufficient to extend this scaling region over a wide 
range of couplings and that as our approximation breaks down we observe a level crossing 
effect. We will examine the possibility of the small (3 minima being scaling regions in a later 
publication. 

The antisymmetric results presented here can also be compared with those of Teper j^J . 
While each of our results is considerably higher than the masses of Teper's and * 
states, Teper's ** state has a mass which is close to the lowest mass state calculated 
here. Teper obtains the following masses, in units of e 2 , for the ** state: 5.42 ± 0.16 
for SU(3), 6.98 ± 0.26 for SU(4) and 9.18 ± 0.45 for SU(5). It is interesting to note that 
our corresponding lowest unimproved glueball masses are consistent with these results. Our 



28 



1.4 H 
1.2 

1 





Unimproved 
Improved 
Tadpole improved 



10 



15 



20 



25 



30 



1.5 h 

1 





Ummproved 
Improved 
Tadpole improved 



10 



15 



20 



25 



30 



(a) Symmetric SU(2) massgap 



(b) Symmetric SU(3) massgap 



t 4. 



5.0 

4.5 

4.0 

3.5 

3.0 

2.5 

2.0 

1.5 




1 1 1 1 1 1 1 






V 




1 1 


Unimproved 

Improved 

Tadpole improved 





.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 

P 



0.5 
6.0 
5.5 
5.0 
4.5 
4.0 
3.5 
3.0 
2.5 
2.0 
1.5 



n 1 1 1 1 1 r 




Unimproved 
Improved 
Tadpole improved 



J I I I I I L 



0.0 10.0 20.0 30.0 40.0 50.0 60.0 70.0 80.0 

P 



(c) Symmetric SU(4) massgap 



(d) Symmetric SU(5) massgap 



FIG. 7: 2+1 dimensional symmetric massgaps for SU(2), SU(3) (both with L r , 
(with L max = 16) and SU(5) (with L max = 12). 



25), SU(4) 



improved results show better agreement with Teper's ** state for SU(3) and SU(5) than 
the corresponding unimproved results. The improved SU(4) results are also consistent with 
Teper's SU(4) ** mass although the agreement is closer for the unimproved result. 

When compared to equivalent unimproved calculations, the improved and tadpole im- 
proved massgaps approach scaling faster as (3 is increased. This is evident in Figs. El and [7| 
and is expected since, for an improved calculation one is closer to the continuum limit when 
working at a given coupling. However, for most improved calculations the scaling behaviour 
is marginally less precise than the equivalent unimproved calculation. A possible reason for 
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TABLE I: Estimates of the lowest lying SU(2) glueball masses (in units of e 2 ) computed with various 
Hamiltonians in 2+1 dimensions. The unimproved, improved and tadpole results are calculated in 
the respective scaling regions 13.5 < < 30.0, 9.9 < < 30.0 and 9.25 < < 30.0. 



Unimproved Improved Tadpole Improved 

Amf 2.05691+0.00002 2.0897+0.0003 2.0965+0.0006 

Amf 3.645+0.001 3.685±0.001 3.6953+0.0009 

Amf 4.5202+0.0004 4.574+0.004 4.583+0.004 

Amf 5.133+0.003 5.177+0.004 5.189+0.004 

Amf 5.867±0.006 5.932+0.008 5.943+0.008 



TABLE II: Estimates of the lowest lying symmetric SU(3) glueball masses (in units of e 2 ) computed 
with various Hamiltonians in 2+1 dimensions. The results are calculated in the scaling region which 
minimises the standard error in each case. 



Unimproved Improved Tadpole Improved 

Amf 3.265868+0.000042 3.32365+0.00012 3.32580+0.00015 

Amf 6.23903+0.00065 6.30391+0.00083 6.31192+0.00084 

Amf 7.5767+0.0025 7.6466+0.0030 7.6498+0.0030 

Amf 8.9462+0.0029 9.0118+0.0044 9.0206+0.0045 

Amf 10.0778+0.0071 10.1546+0.0094 10.1628+0.0094 



TABLE III: Estimates of the lowest lying symmetric SU(4) massgaps (in units of e 2 ) computed 
with various Hamiltonians in 2+1 dimensions. The results are calculated in the scaling region 
which minimises the standard error in each case. 



Unimproved Improved Tadpole Improved 

Amf 4.59121+0.00007 4.6720+0.0001 4.6754+0.0001 

Amf 8.8122+0.0012 8.9276+0.0016 8.9284+0.0017 

Amf 10.5889+0.0051 10.7807+0.0051 10.7794+0.0051 

Amf 12.5527+0.0048 12.6266+0.0081 12.6138+0.0080 

Amf 14.052+0.012 14.165+0.016 14.157+0.016 
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TABLE IV: Estimates of the lowest lying symmetric SU(5) massgaps (in units of e 2 ) computed 
with various Hamiltonians in 2+1 dimensions. The results are calculated in the scaling region 
which minimises the standard error in each case. 



Unimproved Improved Tadpole Improved 

Amf 5.8903±0.0001 5.99434±0.00009 5.9983±0.0002 

Amf Il.2335i0.0036 Il.3696i0.0050 11.3731+0.0049 

Amf 13.340+0.011 13.658+0.012 13.663+0.011 

Amf 15.881+0.012 15.890+0.019 15.890+0.019 

Amf 17.564+0.025 17.676+0.035 17.682+0.035 



TABLE V: Estimates of the lowest lying antisymmetric SU(3) massgaps (in units of e 2 ) computed 
with various Hamiltonians in 2+1 dimensions. The results are calculated in the scaling region 
which minimises the standard error in each case. 



Unimproved Improved Tadpole Improved 

Amf 5.32750+0.00047 5.39661+0.00034 5.39864+0.00032 

Amf 7.9389+0.0021 8.0142+0.0028 8.0145+0.0028 

Amf 8.9319+0.0045 9.0092+0.0056 9.0087+0.0055 

Amf 10.4711+0.0058 10.5514+0.0085 10.5502+0.0085 

Amf 11.304+0.011 11.384+0.015 11.381+0.015 



TABLE VI: Estimates of the lowest lying antisymmetric SU(4) massgaps (in units of e 2 ) computed 
with various Hamiltonians in 2+1 dimensions. The results are calculated in the scaling regions 
which minimise the standard error in each case. 



Unimproved Improved Tadpole Improved 

Amf 7.21479+0.0012 7.3310+0.0011 7.33586+0.00077 

Amf 10.9117+0.0033 11.0099+0.004909 11.0617+0.0046 

Amf 12.121+0.007 12.2779+0.0088 12.292+0.009 

Amf 14.5012+0.0089 14.592+0.014 14.574+0.015 

Amf 15.4521+0.0173 15.545+0.023 15.555+0.023 
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FIG. 8: 2+1 dimensional antisymmetric massgaps for SU(3) (with L max = 25), SU(4) and SU(5) 
(both with L max = 12). 

this is that the one plaquette trial state used here does not allow for direct contributions 
from the improvement term in the kinetic Hamiltonian. For this term to contribute directly 
one would need a trial state which includes Wilson loops extending at least two links in at 
least one direction. 

The improved SU(2) massgap can be compared to the coupled cluster calculation of Li 



et al 36(. Their result (in units of e 2 ), Amf = 1.59, is again significantly lower than our 
result 2.0897 ±0.0003. The difference is again attributable to the different choices of Wilson 
loops used in the simulation of states. While our calculation makes use of the simple one 
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TABLE VII: Estimates of the lowest lying antisymmetric SU(5) massgaps (in units of e 2 ) computed 
with various Hamiltonians in 2+1 dimensions. The results are calculated in the scaling regions 
which minimise the standard error in each case. 



Unimproved Improved Tadpole Improved 

Amf 9.067±0.003 9.239±0.002 9.248±0.002 

Am^ 13.717±0.008 13.89±0.01 13.89±0.01 

Am^ 15.054±0.015 15.32+0.02 15.327+0.019 

Amf 18.08+0.02 18.12±0.03 18.12±0.03 

Am^ 19.084+0.036 19.19±0.05 19.20+0.05 



plaquette ground state and a minimisation basis with only rectangular loops, the coupled 
cluster calculation of Li et al uses a more accurate ground state wave function consisting of 
an exponential of a sum of extended loops which are not necessarily rectangular. Without 
including additional small area Wilson loops we cannot be confident that the lowest mass 
state accessible with our minimisation basis is in fact the lowest mass state of the theory. 
Clearly there is scope for more work here. 

VIII. CONCLUSION 

In this paper we have extended the analytic techniques of 2+1 dimensional Hamiltonian 
LGT, traditionally used for SU(2), to general SU(iV). Impressive scaling is achieved over 
an extremely wide range of couplings for the lowest energy eigenstates in the symmetric 
and antisymmetric sectors. Our calculations use a one plaquette trial state and a basis of 
rectangular states over which excited state energies are minimised. Such choices allow the 
use of analytic techniques in SU(iV) calculations. 

The results of this paper give estimates of the lowest unimproved, improved and tadpole 
improved SU(iV) glueball masses, all of which are above current estimates. We suspect 
that the reason for the discrepancy is a lack of small area non-rectangular states in our 
minimisation basis. A basis of rectangular states was used for simplicity. The inclusion 
of non-rectangular states is straightforward and only complicates the counting of overlaps 
between diagrams of a particular type. When not including sufficient small area diagrams it 
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is possibly that the lowest mass states of the theory are not accessible over a large range of 
couplings. A further improvement to our calculation would involve the use of an improved 
trial vacuum state. Such a state would possibly include several extended Wilson loops in its 
exponent. Without the development of new techniques for performing the required integrals 
the use of such a vacuum wave function would require the use of Monte Carlo techniques for 
the calculation of expectation values. In this scenario many advantages of the Hamiltonian 
approach would be lost. 

In a later publication we extend the calculations presented here to SU(25) in an attempt 
to explore the mass spectrum in the large N limit of pure SU(iV) gauge theory. 
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